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ABSTRACT 

We present sensitive aperture synthesis observations of the nearby, late-type spiral galaxy 
NGC 6503, and produce H I maps of considerably higher quality than previous observations by 
van Moorsel fc Welld (1985). We find that the velocity field, while remarkably regular, contains 
clear evidence for irregularities. The HI is distributed over an area much larger than the optical 
image of the galaxy, with spiral features in the outer parts and localized holes within the HI 
distribution. The absence of absorption towards the nearby quasar 1748+700 yields an upper 
limit of 5 X 10^^ cm^^ for the column density of cold HI gas along a line of sight which should 
intersect the disk at a radius of 29 kpc. This suggests that the radial extent of the HI disk is 
not much larger than that which we trace in HI emission (23 kpc). The observed HI distribution 
is inconsistent with models of a single thin or thick disk. Instead, the data require a model 
containing a thin disk plus a thicker low column-density HI layer that rotates more slowly than 
the thin disk and that extends only to approximately the optical radius. This suggests that the 
presence of extra-planar gas in this galaxy is largely the result of star formation in the disk rather 
than cold gas accretion. Improved techniques for interferometric imaging including multi-scale 
Clean that were used in this work are also described. 



Subject headings: galaxies: individual (NGC 6503) 
— galaxies: spiral — methods: data analysis 

Introduction 



Attention was drawn to the relatively nearby, 
low-luminosity spiral galaxy NGC 6503 by the 
occurrence of a bright quasar only 5.3 arc min- 
utes a way. Deep optical plates taken bv lArp et al. 
(Il976l) show wisps or spiral arms in the outer- 
most portions of the NW part of the galaxy. 



- galaxies: ISM — galaxies: kinematics and dynamics 



and later with the early, incomplete VLA by 
van Moorsel fc Welld ( 19851 ). Neither group found 



any connection between the quasar and the galaxy. 
The images obtained by Ivan Moorsel k WeM 



which were in the direction of, but did not ap- 
proach, the quasar. The neutral hydrogen was 
observed with the W esterb ork Synth esis Ra- 
dio Telescope by IShostak. Willis, fc Crane (,198l|) 
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(|1985I ) revealed a very regular rotation curve, but 
were sensitive only down to a column density of 
1.3 X 10^0 cm-2. 

Although quasars are no longer thought to be 
directly linked to nearby galaxies, NGC 6503 re- 
mains an interesting target: its proximity, isola- 
tion, and regular kinematics make it ideal for de- 
tailed structural studies to understand better the 
formation and evolution of nearby spirals. We 
therefore undertook to re-observe this galaxy using 
the full VLA with improved receivers in the C ar- 
ray in order to obtain substantially better spatial 
and spectral resolution with improved sensitivity. 
Improvements in software allowed us to take ad- 
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vantage of the bright quasar for caHbration with- 
out suffering th e sidelobe problems enco u ntered 
particularly by IShostak. Willis, fc Cran3 (|l98l[ ) 
and further improvements allowed us to image the 
galaxy with very little loss of fidelity due to miss- 
ing short spacings. 

Deep HI observations of nearby spirals have 
revealed a wealth of low column density features 
at anomalous positions and velocities. In par- 
ticular, studies of systems such as NGC 2403 
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(ISchaap. Sanc isi. & Swaters 2000; Fraternali et al. 

20021). NGC 2613 (Chaves fc Irw in 2001), NGC 4559 
dBarbieri et al.l2005h . NGC 89llOosterloo. Fraternah. fc 
20071), U GC 7321 ( ^Matthews fc Wood' '2003") , 
NGC 253 dBosma et al..2005 ). NGC 6946 (B oonisma et al 



bination with high-quality Ha and CO kinematics 
are forthcoming. Here, we describe our HI obser- 
vations and the derived properties of NGC 6503 
in preparation for that work, and model the prop- 
erties of its extra-planar HI. In the appendix, 
we introduce improved image deconvolution algo- 
rithms, describing in some detail the particular im- 
plementation, widely available in the astronomical 
community, used on the present data. 

1.1. Known properties of NGC 6503 



Xhe basic properties of NGC 6503 are sum- 
id in Table [T] The galaxy is a late- type. 



20081 ) and M83 ([Miller et al.l l2008f) suggest that 
extra-planar gas with decreasing rotation veloc- 
ity with increasing height from the disk may 
be ub iquitous in these systems (see iFraternali 
(|2008f) for a complete list). This anomalous 
gas is not clearly associated with previous in- 
teractions or the presence of companions. It 
has been proposed that this extra-planar HI 
stems from a galactic fountain mechanism pro- 
duc ed by feedback from supernovae in the disk 
(e.g'Collins et al."2002l; 'Fraternah fc Binnev"200i 



Oost c rloo, Fratcrna lij fc Sancis i 2007; Boom sma et al 
20081 : iMiller et al.ll2008h . although at least some 



of it must stem from cold gas accretion from the 
intergalactic medium (e.g. .Barnabe et al.. ,2006 



Kaufmann et al. 20061 ; Fraternali fc Binnev 20081 ; 
Sancisi et al.l 120081 ). It has proven difficult to dis- 
tinguish between these two possibilities, however, 
in part because many of the galaxies imaged in HI 
are actively forming stars. 

NGC 6503 is an interesting system to examine 
in this context because it offers the possibility of 
distinguishing between the star formation and gas 
accretion hypotheses. It has a relatively modest 
star formation rate and an HI layer that extends 
far beyond the star- forming disk. If primordial gas 
accretion is a dominant contributor to the extra- 
planar HI, one would expect little correlation be- 
tween this layer and the star-forming disk. On the 
other hand, if star formation were responsible, one 
wouldn't expect to find a thick HI layer beyond the 
optical disk. 

The data set presented in this paper provides 
sensitive HI maps of NGC 6503. Dynamical mod- 
els of this galaxy obtained from these data in com- 



l ow luminosity spira,! locat ed near the Local Void 
( Karachentsev et al. I l2003h . Deep optical surveys 
have not found any galaxies near NGC 6503, mak- 
ing it a very isolated, probably unperturbed sys- 
tem. We note that because of its large angular 
size, it is not included in standard isolated ga laxy 
catalogs (e.g. IVerdes-Montenegro et al.ll2005l ). 

Due to its proximity, the distance to NGC 6503 
cannot be determined reliably from its recessional 
velocity. Photo metry of the brightest s t ars in 
the galax y led iKarache ntsey fc Sharina ( 1997 ) 
also iKarac hcntscv e t al.l ( 20031) ) to derive 



[see 

its distance as 5.2 Mpc, a value that we adopt 
in this study. The Hubble ty pe of SA(s)cd 



from Ide Vaucouleurs et al.l (|l99ll ) indicates that 
NGC 6503 has an S shape and no bar, but 



Lira et al.l (120071) prefer a H ubble type of Sc rather 
than Scd. Makarova (Il999l) presents B, V, and - 
surface brightness profiles traced to 200 arcsec ra- 
dius, and finds the radius of NGC 6503 at 25*'^ 
magnitude per square arcsec to be 116.6, 140.6, 
and 198.9 arcsec, respectively at these bands. The 
central surface brightness of NGC 6503 lies in the 
range found for normal disk galaxies. 

We compute the gl obal star fo r matio n rate of 
NGC 6503 using the iKennicutd (|l998[) IR star 
formation rel a tion using the measurements of 
Sanders et al.l (12003), f o llowin g the method out- 



lined by iKewlev et al.l (|2002n . Its modest star 
formation rate explains its low B-band lumi- 
nosity compared to normal spirals of the same 
morphological type, even when this quantity 
i s normalized by the g alaxy size and/or mass 
(|Roberts fc HavnedlTggl . 

NGC 6503 is in general a fairly normal, if rather 
modest, low- luminosity galaxy with very regular 
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Table 1: NGC 6503 properties 



parameter value reference 



Right ascension (J2000) 


17'H9'"26^30 


this work 


Declination (J2000) 


70°08'40.7" 


this work 


Systemic Velocity (heliocentric) 


28.2 ±0.3 kms-i 


this work 


Distance 


5.2 Mpc 


Karachentsev & Sharina ("1997) 


optical radius (25'^ B mag/ss) 


3.55' 


de Vaucouleurs et al. fl9911 


optical radius 


5.35 kpc 


at 5.2 Mpc 


Hubble type 


SA(s)cd 


de Vaucouleurs et al. (19911 


Absolute B magnitude 


-17.68 


Makarova (1999) 


B luminosity 


1.5 X 10^ Lq 




B central surface brightness 


20.05 mag arcsec"^ 


Makarova ( 1999) 


Star formation rate 


O.18M0yr-i 


this work 



kinematics . Howe ver, in the center of this galaxy, 
Lira et al.l (|2007l ) find a bright central nucleus 
which exhibits "prominent LINER emission lines 
and a red continuum with strong absorption fea- 
tures." Additiona lly, they find six point sources 
of X-ray emission. IHo. Filippenko. SargentI (|l997[ ) 
call the nucleus a "transition" or even Seyfert type 
2 with some uncertainty. The presence of some 
continuum radio emission in a larger, but central 
region of NGC 6503 (Section 13. ip . suggests that 
this galaxy hosts a low-luminosity active galactic 
nucleus (AGN). 

2. Observations and reductions 



NGC 6503 was observed on 27 February 1996 
with the VLA in C configuration. There were 127 
spectral channels separated by 24.4 kHz (5.15 km 
s"""^) centered on 26 km heliocentric radial ve- 
locity. The sources used for bandpass shape and 
absolute flux density calibration were 3C286 and 
3C48, and were observed at the beginning, mid- 
dle, and end of the session for a total of about 43 
minutes. The total integration time on NGC 6503 
was a little more than 500 minutes, interspersed 
every 50 minutes or less with 2-minute observa- 
tions of the 2.0-Jy amplitude and phase calibrator 
1800-1-784. The observing setup is summarized in 
Table H 

All data re duction was d one in the AIPS soft- 
ware package (|Greisenll2003l ). After deleting ^2% 



to NGC 6503 (and 1800-1-784). The bandpass 
shapes were determined from the 3C286 and 3C48 
observations, interpolated in time, and applied to 
NGC 6503 and 1800-f 784. Five channels exhibit- 
ing weak Galactic HI absorption were flagged in 
3C48, and the bandpass was interpolated across 
them before being applied to the data. 

The 0.76-Jy quasar 1748-1-700, in the same 
beam as NGC 6503, dominates the data in all con- 
tinuum channels. To separate the continuum and 
line signals, the data for NGC 6503 were phase 
shifted to the quasar position and straight lines 
(with slope) were fit to the spectrum of the real 
and imaginary parts of each visibility sample us- 
ing a total of 50 channels from the band edges 
which were judged to be both free of line signal 
and bandpass edge effects. A line-only dataset was 
produced by subtracting the fits from the visibil- 
ity data and shifting back to the original phase 
center, while the best-fitting continuum values at 
the band center were adopted as the continuum- 
only dataset . A task similar to the on e used is 
described bv lCornwell. Uson. &: Haddad (|l992h . 

All imaging and deconvolution were done with 
the AIPS program IMAGR using a t empered uni- 



form w eighting ( "robust weighting" , iBriggs et al 
jiggi). The point-source response was nearly cir- 



of the data for various reasons such as scan start- 
up problems and spurious amplitudes and phases, 
the amplitude and phase calibrations were applied 



cular. We used IMAGR's multi-scale Clean algo- 
rithm to deconvolve the synthesized beam from 
both the continuum and line datasets. Because 
this functionality has not yet been thoroughly de- 
scribed in the literature, we discuss it in detail in 
Appendix El The final image properties are sum- 
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marized in Table [21 

The continuum was imaged over an area ap- 
proximately 2.5° in diameter in order to cover 
both the primary beam and the first outer side- 
lobe of the individual VLA antennas. A few more 
distant regions surro unding nearby NVSS sources 
( Condon et al. |[l99i) were also included. The area 
was imaged in 55 facets (|Cornwell fc Perlevlll992r ) 
in order to avoid wide-field distortions. The decon- 
volution was performed with a multi-scale Clean 
using Gaussian model sources with full-width half- 
maximum (FWHM) sizes of 0", 36", and 108", and 
Cleaning was stopped when the peak residual in 
each facet fell below a flux density cutoff of 0.21, 
0.42, and 1.7 mJy/beam, respectively. The "bias" 
parameter controlling the facet Cleaning order was 
b = 0.62 (see Appendix [X]). 

The proximity to NGC 6503 of the bright 
quasar 1748-1-700, as well as other continuum 
sources, allows for an accurate self-calibration to 
improve the amplitude and phase ca, libration, as 
initially described by ISchwabI ( 1980l ). The im- 
plementation in AIPS uses the full multi-scale, 
multi-facet Clean model to compute model vis- 
ibilities which are divided into observed visibili- 
ties to get apparent baseline-dependent gains. A 
least-squares method turns these into antenna- 
dependent gains which are applied to the observa- 
tions before another round of imaging. Given the 
quality of the initial calibration using 1800-1-784 
and the strength of 1748+700, the self-calibration 
iterations converged quite rapidly. After imag- 
ing, deconvolution and several rounds of self- 
calibration, the 55 separate facet images were 
interpolated onto a "flattened" image with no 
loss of intensity or positional accuracy. The fi- 
nal continuum image has a resolution of 14", a 
root mean square (RMS) of 80 /iJy/beam and a 
dynamic range in excess of 9000:1. 

The improved amplitude and phase solutions 
from the self-calibration and some additional edit- 
ing were then applied to the HI visibility data. 
Each channel was imaged in a single 50' x 50' 
facet. The deconvolution was performed with a 
multi-scale Clean using Gaussian model sources 
with FWHMs of 0", 36", 108" and 324", flux den- 
sity cutoffs of 0.2, 0.55, 2.6, and 7.8 mJy/beam 
respectively, and a bias parameter b = 0.62 (see 
Appendix[X]). The Clean was limited interactively 
to the regions which clearly had emission in that 



chan nel in order to mitig ate the well-known Clean 
bias (jCondon et al.lll998[ ). The final HI data cube 
has an RMS noise of 0.58 mJy/beam, a spatial 
resolution of 14" per channel, and a peak dynamic 
range of 30:1. 

3. Images 

The properties of the full-resolution continuum 
and HI images of NGC 6503 are summarized in 
Table [21 Below, we discuss each in turn. 

3.1. Continuum in NGC 6503 and in 
1748+700 

The total Cleaned flux density from the 2.5° 
diameter continuum image in the vicinity of 
NGC 6503 is 967 mJy, uncorrected for attenu- 
ation from the primary beam of the individual 
VLA antennas. There are ~ 90 sources with a 
signal-to-noise ratio S/N > 5 in the image, some 
of which are in the first and second outer sidelobes 
of the primary beam. The sources are essentially 
all unresolved except for NGC 6503. 

Fig- dl presents a portion of the primary beam- 
corrected image of the continuum emission in the 
vicinity of NGC 6503. As expected, the noise in 
the image rises rapidly towards the edge of the 
primary beam. The total image size is governed 
by the beam correction itself, since the VLA's pri- 
mary beam is completely unknown beyond 30' at 
this frequency. The paucity of continuum sources 
near NGC 6503 and 1748+700 as well as the ap- 
parent chain of sources east of the pair is initially 
striking, but not statistically significant (J. J. Con- 
don, priv. comm.). The extent of the HI emission 
in NGC 6503 (see Section 14. ip is represented by 
the solid line in Fig. [l] Qualitatively, the centroids 
and position angle of the HI and continuum emis- 
sion are in good agreement, with the HI extending 
well beyond the continuum. 

The continuum properties of NGC 6503 are 
summarized in Table [3l The primary-beam cor- 
rected total flux density of the extended contin- 
uum source in NGC 6503 is 39 ± 1.2 mJy. A 
contour image of the continuum radio source in 
NGC 6503 is shown in Fig. [51 The position angle 
is close to the -60 degrees found in the HI. The full 
extent along the major axis is about 200"or 5 kpc 
at the adopted distance. The continuum appears 
to peak at the three locations marked by crosses. 
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Fig. 1. — Primary beam-corrected continumn im- 
age of a field near NGC 6503. The wedge to the 
right shows the displayed brightness in mJy/beam 
which is clipped at 5 mJy/beam. The extent of the 
detected HI emission is shown by the solid line (see 
Section HH]). The quasar 1748-1-700 is the strong 
point source 5.34' to the SW of the galaxy. 



Table 2: Observing setup and image properties 



parameter 


value 


units 


VLA array configuration 


C 




Primary beam FWHM 


29.7 


arc minutes 


Time on-source 


500 


min. 


Usable total bandwidth 


2513 


kHz 


Band center (heliocentric) 


26 


km 


Spatial resolution 


14 


arcsec 


Spatial resolution 


0.35 


kpc 


Effective cont. bandwidth 


1220 


kHz 


RMS image noise (cont.) 


0.08 


mJy/beam 


Dynamic range (cont.) 


9000 




Channel separation (line) 


5.15 


km s~^ 


RMS image noise (line) 


0.58 


mJy/beam 


Peak dynamic range (line) 


30 





o 




RIGHT ASCENSION (J2000) 



Fig. 2. — 21-cm radio continuum distribution in 
NGC 6503. The wedge at the top shows the dis- 
played brightness in mJy/beam, and the contours 
are drawn at (1, 2, 3, 4, 5) x 0.25 mJy/beam. The 
Clean beam is plotted in the top left corner. The 
continuum peaks are marked by crosses, th e 6 un- 
resolved X-ray sources found bv iLira et al. ( 2007 ) 
are marked by stars, and the kinematic center of 
the HI distribution (Section 14. 2p is marked by a 
plus sign. 
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The HI is at a minimum in the center, but has 
peaks near the other three locations. The appar- 
ent ridge in the continuum in the NW crosses a 
region in which the HI is minimal, suggesting that 
this correspondence is probably not c ausal. Th e 



six un resolved X-ray sources found by iLira et al 
POOTI ) are also illustrated with stars. That the 



central continuum radio peak is coincident with 
an X-ray source is further evidence for a low - 
luminosity AGN (jHo. Fihppenko. Sargeii^ll997l) . 

With primary beam correction, the total flux 
density of the quasar 1748-f 700 was 755 mJy on 
the observation date. Within reasonable uncer- 
tainties the quasar is unresolved. 

3.2. HI in NGC 6503 

Fig. [3] shows the detected HI emission from 
NGC 6503 in a selection of channels. It is clear 
that even at the higher resolution and sensitiv ity 
of our data relative to those of lvan Moorsel fc Wellsl 
(|l985l) . the HI kinematics of NGC 6503 are re- 
markably regular. Nonetheless, when the channel 
images are examined in detail, irregularities are 
apparent. The outer portions of many emission 
"wings" bend away from the major axis, most no- 
ticeably at -25.5 km s~^, -10.1 km s~^ and 82.7 km 
in Fig. [HI There are also spatially distinct fea- 
tures in the distribution of some channels, such as 
at -25.5 km s~^ and 98.2 km s~^. We investigate 
the HI morphology and kinematics of NGC 6503 
in detail in Section [H 

Fig. [4] shows selections from the data cube af- 
ter rotation and transposition, plotting the HI 
as a function of velocity and position along lines 
parallel to the major axis. The central panel 
is the traditional rotation curve image of this 
galaxy. Fig. [5] shows position-velocity images 
along the minor axis at several positions on the 
major axis. These figures are made from the full 
resolution image using a single 3-arcsecond pixel 
on the third axis for each panel with no addi- 
tional smoothing. At first glance, the position- 
velocity slices of NGC 6503 exhibit considerable 
symmetry and do not appear to show an exten- 
sive "beard" of HI emission at velocities between 
the rotation velocity and the syst emic velocity, as 
found for several other systems dFraternali et al 



2002"; 'Matthews & Wood''2003'; 'Bosnia eLalJ|2005|; 

Oosterloo, Fratcrnali, & Sancisi 2007). In Sec- 
tion [431 we place quantitative constraints on the 



morphology and kinematics of the extra-planar HI 
layer that is consistent with the present data. The 
bottom panel in Fig. [5| in particular shows some 
significant asymmetry. This is associated with the 
warp to the SW discussed further in Section [¥751 in 
connection with Fig. [131 

3.3. HI in the vicinity of NGC 6503 

We searched the data cube outside the confines 
of NGC 6503 for any other source of line emission. 
Due to limitations of the VLA correlator and the 
dishes that make up the array, a bandwidth of only 
about 500 km s~^ was observed, centered at 26 km 
s~^ heliocentric, over a region about 50'in diame- 
ter. No other extragalactic HI emission was found. 
For an unresolved source spanning two contiguous 
velocity channels at the distance of NGC 6503, this 
non-detection corresponds to a 3(7 HI mass limit 
that varies between Mhi = 1.1 x 10^ Mq at the 
pointing center and ten times this value near the 
edge of the field. In the limited volume probed by 
our observations, NGC 6503 appears isolated from 
an HI standpoint as well as from an optical one. 

3.4. Absorption towrards 1800+784 and 
1748+700 

The calibration source 1800+784 is located at 
Galactic longitude / = 110.0°, latitude b = 29.1°, 
while NGC 6503 is at Z = 100.5°, b = 30.7°. These 
directions are expected to have local Milky Way 
HI emission and, most likely, absorption. Indeed, 
Fig. [6l shows an absorption line in the spectrum of 
1800+784, with a depth of 0.040 ± 0.002 Jy or an 
optical depth of 0.02. The line is centered at an 
LSR velocity of -1.6 km s^^ with a FWHM of 13 
km s""'^. 

By contrast, Fig. [7l shows that no absorption is 
observed in the direction of the quasar 1748+700 
{I = 100.5°, b = 30.8°) to very high precision. 
The deepest negatives in the full resolution im- 
age in this direction (light line) are less than 1 
mjy/beam, while the spectrum extracted from an 
image smoothed to 30" and 10 km resolution 
(heavy line) shows even less evidence for absorp- 
tion. The optical depth toward 1748+700 is thus, 
conservatively, r = 0.0+0.001. This suggests that, 
for HI of line width 5 km s~^ and temperature 50 
K, the column density upper limit is 5 x lO^^cm"^. 
Assuming a thin, flat disk in NGC 6503, the line 
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Fig. 3. — Primary-beam corrected channel images of the HI emission in NGC 6503. The gray-scale is plotted 
from 0.5 to 18.65 mJy/beam with a logarithmic transfer function and contours are plotted at -2, 2, 4, 8, and 
16 mjy/bcam. (The negative contours appear only in regions which arc noise dominated.) The heliocentric 
radial velocity of the channel is in the upper right corner, and the Clean beam is plotted in the lower left 
corner of the first panel. The kinematic center of NGC 6503 is marked with a plus sign and the quasar 
1748-1-700 is marked with an asterisk. 
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Fig. 3. — Continued. 
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Fig. 4. — Primary-beam corrected images for slices parallel to the major axis, separated by 21". Axis labels 
are relative to the kinematic center and the systemic velocity. The gray-scale is plotted from 0.5 to 18.65 
mJy/beam with a logarithmic transfer function and contours are plotted at 1.5, 3, 6, and 12 mJy/beam. 
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Fig. 5. — Primary-beam corrected images for slices parallel to the minor axis, separated by 42". Axis labels 
are relative to the kinematic center and the systemic velocity. The gray-scale is plotted from 0.5 to 18.65 
mJy/beam with a logarithmic transfer function and contours are plotted at 1.5, 3, 6, and 12 mJy/beam. 
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Fig. 6. — The spectrum of calibration source 
1800+784 computed by vector averaging all cal- 
ibrated fringe visibilities. 
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Fig. 7. — The spectrum at the position 
of 1748+700 (RA = 17'^48'"32!84, Dec = 
+70°05'50'.'8 in J2000) minus the measured con- 
tinuum flux density of 755 mJy. Light line from 
the high-resolution, primary-beam corrected data 
cube (RMS = 0.58 mJy/beam). Heavy line from 
a version of the data cube smoothed to 30" spa- 
tial resolution and 10 km s~^ velocity resolution 
(RMS = 0.38 mJy/beam). 



of sight to the quasar passes through the plane of 
NGC 6503 at a radius of 29 kpc. There must be 
very httle cold HI at this radius in NGC 6503. The 
absence of absorption from the Milky Way is also 
remarkable. 

We note that in Fig. [71 the peak at an LSR 

velocity of 30 km s"! with a FWHM ~ 40 km 

s~^ coincides with the channels contaminated by 
Milky Way emission over much of the field. The 
projected heliocentric velocity of the HI disk in 
NGC 6503 along this line of sight is ~ 50 km s~^ 
(LSR velocity of ~ 65 km s~^); we therefore do 
not expect Milky Way emission or absorption to 
confound an absorption signature from NGC 6503. 

4. HI in NGC 6503: analysis 

In this section we examine the HI morphology 
and kinematics of NGC 6503 in detail. 

Analysis of an image over an area much larger 
than that containing real signal produces large 
uncertainties due to the inclusion of regions con- 
taining only noise. To mitigate this effect in 
NGC 6503, we produced a cube in which line-free 
regions were blanked. The multi-scale Cleaned im- 
age cube was smoothed in frequency with a Gaus- 
sian of 3 channels FWHM and convolved to a spa- 
tial resolution of 25", producing an image with 
RMS 0.22 mJy/beam. Then all pixels in the fuU 
resolution cube for which the smoothed cube was 
less in absolute value than three times the RMS of 
the smoothed cube were blanked. An interactive 
display program was used to eliminate the remain- 
ing pixels throughout the cube that were discon- 
nected from each other and from any region of 
obvious line signal. 

4.1. Total HI and its distribution 

The blanked image cube may be analyzed along 
the velocity axis for its zeroth, first, and second 
moments. The zeroth is a measure of total HI con- 
tent while the first, if the spectra are narrow and 
symmetric as they are in NGC 6503, is a measure 
of the local velocity. Under the same conditions, 
the second moment reflects the velocity dispersion 
of the gas. The first two moments are shown in 
Fig.m while the second moment is shown in Fig. [9] 
In this section, we focus on the features of the ze- 
roth moment map. The HI properties derived from 
this distribution are summarized in Table [3l 
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Table 3: Measured properties of NGC 6503 



parameter 


value 


units 


Section 


Continuum flux density 


39.4 ±0.7 


mJy 


[Q 


Continuum angular diameter 


200 ± 20 


arcsec 


[Q 


Continuum linear diameter 


5 ±0.5 


kpc 


[Q 


Line integral 


205 ± 1 


Jy km s~^ 


1411 


HI mass 


1.3 ±0.2 


IO^Mq 




HI mass / B luminosity 


0.87±0.13 


Mq/Lq 




Integrated line width at 20% peak 


246 


km 




Integrated line width at 50% peak 


234 


km 




HI radius (0.1 Mq /pc^) NW 


22.6 


kpc 




HI radius (O.lMo/pc^) SE 


17.2 


kpc 


KT\ 


Mean Hi/optical radius ratio 


3.7 






Kinematic center (RA, J2000) 


17h49"*26".30 




14:21 


Kinematic center (Dec, J2000) 


70°08'40.7" 




14:21 


Kinematic Position angle 


-60.1 


degrees 


14:21 


Inclination 


75.1 


degrees 


01 


Systemic velocity (heliocentric) 


28.5 ±0.3 


km 


K2\ 


dynamical mass {R = i?max= 3.5 kpc) 


1.0 


10^° Mq 




dynamical mass {R — 6 kpc) 


1.8 


lOi" Mq 


10 


dynamical mass {R = 20 kpc) 


6.0 


10^° Mq 


10 


Max. column density** 


6.9 


Mq/pc^ 


01 


Max. column density** 


8.6 


10^" atoms cm^^ 


1431 


Minimum detectable signal 


3 


Jy/beam m s^-'^ 




Minimum detectable column density 


1.8 X IQis 


cm^^ 





"corrected for inclination (75.1°) 
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Fig. 8. — Top: zeroth moment of the HI distribution in NGC 6503 measured in 10^^ cm~^ and plotted with 

a logarithmic gray-scale. The column density is 0.57 x 10^^ cm^^ SdV , with SdV in units of Jy/beam m 
s~^.) The lowest reliable HI emission from NGC 6503 is at 1.8 x 10^^ cm^^. Bottom: first moment of the 
HI distribution in NGC 6503. The heavy contour is plotted at 26 km heliocentric (nearly the systemic 
velocity), and the other contours are spaced at 20 km s~"^ relative to this value. The eastern (E) side is 
approaching, while the western (W) side is receding. 
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Fig. 9. — Second moment of the HI distribution in 
NGC 6503 plotted with a hnear gray-scale from 
to 39 km s~^. No correction for projection effects 
has been made. 
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Fig. 10. — The integrated spectrum of NGC 6503 
after primary beam correction and blanking of 
line-free regions. The LSR as well as heliocen- 
tric velocities are shown, with a long tick mark at 
km s~^ LSR. Errors in the spectrum due to HI 
emission from the Milky Way are evident near this 
velocity. 



The weakest pixels in the zeroth moment image 
in Fig. [5] have values as low as 1 Jy/beam m s""'^, 
but the spectra at these pixels do not contain con- 
vincingly real line signals even in the smoothed im- 
age cube. Convincing line signals are seen at pixels 
with zeroth moments of 3 Jy/beam m s~^, yielding 
the lowest reliable HI detection from NGC 6503 of 
1.8 X lO^^ cm-2. 

Integrating over the zeroth moment image 
yields an integrated flux 5int = 205 ± 1 Jy km 
s~^. The corresponding HI mass is given by 

A/hi = (2.356 x 10^)SintD^ Mq , (1) 

where S'int is in Jy km s^^ and D is the dis- 
tance to NGC 6503 in Mpc. The total HI mass of 
NGC 6503 is therefore Mm = 1.3 ± 0.2 x 10^ M©, 
including a 10% uncertainty in D. The integrated 
spectrum of NGC 6503, obtained by summing the 
flux in each channel of the blanked image cube, 
is shown in Fig. [TO] and yields the same HI mass 
within the uncertainties. The width of the inte- 
grated profile is Wsini = 234 km at the half- 
peak points (246 km s~^ at 20% points). 

Several features should be noted in the zeroth 
moment image in Fig. [5] The first is the appear- 
ance of spiral structure, particularly in the outer 
parts of the disk. The approaching arm along 
the southeastern (SE) edge is particularly promi- 
nent. Other arms are visible in the interior to 
the SE, NE, and NW. The spiral arms also ap- 
pear conspicuously in residual images after fits to 
a smooth HI distribution by CUBIT (Section EM 
and Fig. [TO)) . Note that, if the spiral arms are 
trailing arms in the usual sense, the southern edge 
of NGC 6503 is closer than the northern edge. 

There are also a number of holes in the HI dis- 
tribution in Fig. [51 where a distinct minimum is 
surrounded by a ridge of HI. One of these is at 
the very center of the NGC 6503 with an appar- 
ent FWHM of 20". Deconvolved from the beam, 
this diameter is about 14" or 0.35 kpc. This hole 
is visible in the channel and transposed images of 
Figs. [3| through [5] Presumably either the neu- 
tral hydrogen has been converted to molecular 
form or it has been ionized or evacuated by the 
activity at the center of the galaxy which gives 
rise to the weak continuum radio source. South 
of the center there is another clear depression 
with an apparent width of 17", or a deconvolved 
width around 0.24 kpc. Such HI bubbles are seen 
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in other nearby galaxies: for example, detailed 
observ ations of M33 reveal a multi t ude o f these 
holes ( Thilker. Braun. fc WalterbosI liooi) which 
appear to be correlated wi th large OB associations 
(|Deul fc den HartodllQiol ). 



Shostak. WilHs. fc Cran3 (|l98ll ) found that 



NGC 6503 is considerably more extended to the 
SE than to the NW, with 1.25 times more HI on 
the E or approaching half. Measuring the zeroth 
moment image in some detail, we do not confirm 
these early results. Instead we find that the aver- 
age HI in the approaching half is about 11% more 
per pixel but over 8% fewer pixels. The net is only 
4% more HI in the approaching side of the galaxy. 
The outer major axis radius to the NW is 900" 
(22.6 kpc), while it is only 825" (17.2 kpc) to the 
SE. 

The mean Hi/optical radius ratio is therefore 
3.7 for NGC 6503. While this value is larger 
than typically observed in spirals teroeils fc Rhee 
1997t ). it does not rival the relative HI extents of 
some low-mass systems (e.g. Gentile et all ( 20071 )). 
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Fig. 11. — Rotation curve for the HI disk of 
NGC 6503, obtained from the MPS program GAL 
by fitting all points in the first moment image. 



4.2. Kinematics 



The first moment of the HI distribution in 
Fig. [8] shows that the kinematics of NGC 6503 
are highly regular. We find some evidence for a 
mild warp in the outer disk, fi rst identified by 
Shostak. Willis, fc Grand (|l98l[ ): the isovelocity 
contours in Fig. [8] make a somewhat greater angle 
with the major axis in the NE and SW quadrants 
compared to the NW and SE ones. The emission 
in some individual channels of the image cube ex- 
hibits the same behavior, notably at -10.1 km s~^ 
and 67.2 km s^^ in Fig. [3] 

The second moment image (Fig. |9]) has larger 
values toward the center of the galaxy where the 
steepness of the rotation curve produces an in- 
crease in the line width due to the finite beam 
width. There is also a quite noticeable peak in 
the second moment near where the SE arm bends 
abruptly. That feature is due to the presence of a 
second spectral feature in that area, slightly offset 
in velocity from the main feature; see the emission 
at -54.4 km s~^, -41.0 km s~^ and -25.5 km s~^ in 
Fig. [3l It is plausible that the bend in the isoveloc- 
ity contours in the first moment at that location 
(Fig. [HI) also stems from this limited spectral fea- 
ture. 
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Fig. 12. — The rotation curve fit by the GAL 
program for the NW (cross symbols), SW (trian- 
gles), NE (circles), and SE (squares) quadrants. 
Note that the velocity tick marks are separated 
by about one spectral channel in this blown-up 
display. 
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Van Moorsel & Wells (1985) describe an AIPS 
program GAL which fits rotation curves to user- 
selected portions of the first moment image with 
optional weighting by the zeroth moment image. 
We first used GAL to model the kinematics of the 
NGC 6503 HI disk as a whole: after fitting for the 
kinematic center of the disk, this value was held 
fixed and the position angle of the receding major 
axis and the systemic velocity were derived. The 
best fitting values for these parameters are given 
in Table [3] Holding these parameters fixed, we 
then proceeded to fit for the inclination and the 
rotation curve of the system using the functional 
form: 



V = Vm 



-ln(100)i?,/fl„ 



(2) 



where Knax is the maximum velocity of the rota- 
tion curve, which occurs at a radius i?max- 

The rotation curve obtained by fitting all points 
in the first moment image is given in Fig. Illl and 
the extracted parameters are given in the first row 
of Tabled) The rotation curve in Fig. [TT] rises out 
to i? = 140" (3.5 kpc), and then remains flat out 
to i? = 770" (19.4 kpc). There is some evidence 
for a decrease in rotation amplitude beyond this, 
out to the last measured point at i? = 900" (22.6 
kpc). Given a rotation curve, the dynamical mass 
enclosed within a radius R can be computed as- 
suming spherical symmetry: 



M^ 



(2.325 X IQ^)RVL Mr.- 



(3) 



where R is in kpc, and Vmt is the rotation curve 
amplitude in km s^^ measured at R. Dynamical 
mass estimates at different R for the rotation curve 
in Fig. [11] are given in Table [3l 

We then extracted the rotation curves and incli- 
nations of NGC 6503 in the NW, NE, SW and SE 



Table 4: GAL fit parameters by quadrant 



quadrant 


inclination 




RmaK 




deg 


km s^^ 


arcsec 


all 


72.6 


113.2 


139.2 


SW 


71.7 


111.4 


129.7 


NW 


74.3 


111.3 


135.5 


NE 


72.7 


114.5 


149.9 


SE 


72.3 


117.1 


145.7 



quadrants of the disk defined by the kinematic ma- 
jor and minor axes. The results from these fits are 
shown in Fig. [12] and Table [4] Despite the regular- 
ity of the first moment image for NGC 6503, the 
discrepancies between the best fitting parameters 
for each quadrant suggests the presence of non- 
circular motions and geometric distortions in the 
disk. The SW quadrant inclination differs from 
the NW by 2.6 degrees in the direction of the plane 
bending down in that quadrant. The eastern side, 
however, gives an intermediate inclination which 
is the same in both quadrants, failing to confirm 
the second half of the warp suggested above. We 
note that the velocities and physical dimensions 
in the eastern half are somewhat larger than in 
the western half. Fig. [12] makes it clear that there 
are small-scale deviations, of order plus/minus one 
spectral channel (5 km s~^), from simple rotation 
in all four quadrants. The deviations are great- 
est in the SE and NW in which the view is along 
spiral arms against the direction of motion in the 
SE and in the direction of motion in the NW. Pe- 
culiar velocities due to spiral structure would be 
expected in those locations. 

4.3. CUBIT models 

Another AIPS task CUBIT, was initially de- 
veloped by IrwinI (1994) to fit the entire data cube 
rather than just the first moment image. The 
model assumes that the galaxy is a disk of finite 
thickness in circular rotation, and fits for one of 
several functional forms of the HI density distri- 
bution and rotation curve. The CUBIT model fit 
to NGC 6503 held the kinematic center and sys- 
temic velocity fixed to the values in Table [3] while 
fitting for a Gaussian density profile within the 
plane of the disk, an exponential density profile 
perpendicular to the plane, and a Brandt rotation 
curve given by 



V{R) 



Vb R/Rb 



1/3 + 2/3(i?/i?B)" ) 



3/(2m) 



(4) 



where Vb and Rg characterize the rotation curve 
shape and m is the "Brandt index." The results of 
the CUBIT fits to the NGC 6503 cube are summa- 
rized in Table [5] where the columns labeled "mini- 
mum" and "maximum" illustrate the range in pa- 
rameter values found when CUBIT is restricted to 
different regions of the disk. In general, there is 
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Fig. 13. — Top: zeroth moment image of NGC 6503 minus that fit by CUBIT plotted with a logarithmic 
gray-scale. The step wedge is marked in steps of 10^^ cm^^. A comparison of the CUBIT model and the 
radial column density along the major axis is shown in Fig. 1141 Bottom: first moment image of NGC 6503 
minus that fit by CUBIT. The white contour is drawn at km s^^, and the black, dashed contour is drawn 
at -20 km s"^ 
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Table 5: Fit parameters from CUBIT 



parameter 


units 


lull galaxy 


minimum 


maximum 


Position angle 


degrees 


-60.3 


-62.4 


-58.2 


Inclination 


degrees 


75.1 


73.5 


76.0 


Vb 


km s~^ 


118.7 


118.3 


119.9 


Rb 


arc sec 


172.6 


168.5 


218.7 


Brandt index 




0.091 


0.060 


0.69 


Maximum density 


cm~'^ 


0.147 


0.134 


0.180 


Radial density scale*^ 


arcsec 


353.5 


329.2 


373.3 


Vertical density scale*^ 


arcsec 


26.0 


23.6 


33.0 


Vertical column density 


1020cm-2 


4.1 


4.1 


4.8 



"the radius or height at which the density falls by a factor of 1/e = 0.3679. 



good agreement between the best fitting parame- 
ters found by GAL (Table H and by CUBIT (Ta- 
ble E]) . We note that a higher inclination is ex- 
pected from CUBIT, since GAL assumes a thin 
disk while CUBIT takes its finite thickness into 
account. 

Fig. [13] shows the residuals when the best fit- 
ting CUBIT model for the disk as a whole is sub- 
tracted from the data. The residuals emphasize 
the geometric and kinematic deviations from an 
ideal disk that were discussed in Sections 14.11 and 
14.21 They are largest along the SE major axis, 
where the CUBIT model greatly underestimates 
the zeroth moment and does not fit the velocity 
profile well. The kinematic discrepancy between 
the data and the model might be expected given 
the rotation curve fits in Fig. [121 where the largest 
disparity between the quadrants within the opti- 
cal radius is between the SE and SW at i? ~ 130". 
The first moment velocities in the data also exceed 
the model by up to 20 km s~^ along the southern 
edge. This provides some support for the previous 
suggestion of a warp in the SW. 

Integrating the exponential vertical density 
profile, the CUBIT model provides the best fit 
Gaussian to the radial column density profile in 
NGC 6503. This profile can also be estimated di- 
rectly from the zeroth moment image. Different 
inclination-corrected measures of the radial col- 
umn density in NGC 6503 are shown in Fig. [TU 
where the plus signs show the CUBIT model out- 
put, the squares show the profile obtained by av- 
eraging concentric rings in the zeroth moment 




Kilo pc in radius 



Fig. 14. — The radial column density in Mo/pc^ 
plotted against radius in NGC 6503. The plus 
signs trace the CUBIT model, the squares are av- 
erages over concentric rings in the zeroth moment 
image, and the dashed and solid lines are slices 
in the zeroth moment image along the major axis 
to the SE and NW, respectively. The slices have 
been corrected to radial column densities by mul- 
tiplying by the cosine of the inclination (75.1°). 
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image, and the dashed and solid lines show the 
profiles obtained from slices of the zeroth moment 
image along the major axis to the SE and NW, 
respectively. A comparison between the CUBIT 
model output and the SE and NW major axis 
profiles explains the large zeroth moment residu- 
als in Fig. [131 the surface density profile is well- 
approximated by a Gaussian along the NW major 
axis, but deviates from this functional form by 
almost a factor of two along the SE major axis 
at i? '--^ 2kpc (80"). As a result, the zeroth mo- 
ment residuals in Fig. [13] contain nearly 45% of 
the total HI in NGC 6503. We note that the radial 
column density profile estimated from the zeroth 
moment image is contaminated by the relatively 
poor resolution along the minor axis due the high 
inclination of NGC 6503: this biases the profile 
low at intermediate R and high at large R. 

The peak column densities and other param- 
eters found here ar e similar to those found by 
van der Hulst et al. 1 ([1993) 

for a number of low 
surface brightness galaxies. Like them, NGC 6503 
falls below the average for field Scd galaxies (8.5 
MfT) / p c^) and would fall below the iKe nnicuttl 
()l989l) criterion for normal star formation. 



data, we smoothed the HI cubes presented in Sec- 
tion [321 spatially and spectrally to a resolution of 
30" and 10.3km s^^. We compare our models to 
this smoothed cube. 



4.4. Extra-planar HI Models 

In comparison to th e HI distributi ons of 

several nearby spirals ( Fraternah et al.l 20021: 



We use the GIPSY (|van der Hulst et al.lll992D 
task GALMOD to model NGC 6503 as one or 
two HI layers of finite thickness in circular rota- 
tion, using the results from the analysis in Sec- 
tions 14.21 and 14.31 as inputs. In all of the models, 
we adopt a kinematic center, position angle and 
systemic velocity identical to the values in Table [3l 
and a constant velocity dispersion of 8 km 
(Fig. [9]). 

Our CUBIT models suggest that the receding 
(W) side of the disk is better approximated by 
a smooth HI layer than the approaching (E) side 
(Fig. [13]) , and we therefore focus on reproducing 
its properties with GALMOD. Specifically, the 
input model kinematics are the rotation curve ex- 
tracted by GAL for the receding half of the disk, 
and the input HI surface mass density is that along 
the NW major axis (Fig. [T4]l . We have corrected 
the input rotation curve and HI surface density for 
beam smearing at i? < 100". The output models 
are smoothed to the same spatial and spectral res- 
olution as the smoothed data cube. 

A selection of galaxy models is compared to 
the smoothed NGC 6503 image cube in Fig. [15] 
The left column shows a position-velocity slice 



Matthews fc Woodl[2003l:lBosma et al.ll2005l:lOosterloo. FrgfttenmlthfeSfficiBlijor axis, and the right column 



20071) . position-velocity slices along the major 
(Fig. H and minor (Fig. [5]) axes of NGC 6503 
show little evidence for low-velocity gas and a 
high degree of symmetry, respectively. Our CU- 
BIT models of NGC 6503 (Section US]) suggest 
that the mean vertical density scale0 of its HI layer 
is 26" (0.65 kpc). Taken together, the position- 
velocity slices and CUBIT models thus suggest 
that NGC 6503 lacks the thick, lagging extra- 
planar HI layer that is commonly detected in spiral 
galaxy observations at this sensitivity (see Section 
1 for references). In this section, we model the 
vertical extent and kinematics of the HI layer in 
NGC 6503 to place quantitative constraints on the 
properties of its extra-planar gas. 

To emphasize any large-scale, low column- 
density features that might be present in the 



^height over which the density drops by a factor of 1/e 
0.3679 



shows the detected HI emission in the channel at 
V ~ 103.3km s^^. The contours in all rows of 
a column are identical. We choose to show these 
model properties because we find that they pro- 
vide the best diagnostics of the vertical extent and 
kinematics of the HI layer. We note that because 
of projection effects there was little difference be- 
tween the position-velocity slices along the minor 
axis for various models, and they are not included 
in Fig. 

We first model the HI in NGC 6503 as a single 
Gaussian layer. Figs. [TSb and[l5b show the model 
morphologies for HI layers with FWHM = 0.5 kpc 
and FWHM — 4.5 kpc, respectively. They illus- 
trate that a model with a thin HI layer (Fig. [I5b) 
reproduces the gross features of the major axis 
position-velocity slice beyond the optical radius 
(i? = 3.55') as well as the observed channel maps: 
this is expected given the low vertical density scale 
returned by CUBIT (Section [43]). However, this 
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Fig. 15. — Position-velocity slices along the NW major axis (left) and channel map at = 103.3km s"-*^ 
(right) in a) the smoothed image cube, and b) - e) various GALMOD models. The contours in the position- 
velocity panels are at (1, 3, 6 , 10 , 18 , 28 , 40) x 0.76 mJy/beam, and the contours in the channel map panels 
are at 1.5, 4 ,8 ,16 and 24 mJy/beam. Panels b) - e) show our most plausible models including b) a single 
thin HI layer, c) a single thick HI layer, d) thin and thick layers with a column density contrast of 4:1 where 
the thick disk lags the thin one by 25 km s~^, e) thin and thick layers with a column density contrast of 4:1 
where the thick disk lags the thin one by 25 km s~^ and is truncated at the optical radius {R = 3.55'). 
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model is unable to reproduce the position- velocity 
contours within the optical disk: there is more 
anomalously low-velocity gas at these radii than 
expected. Fig. [T5b demonstrates that invoking a 
thick HI layer does produce low-velocity gas at 
small R in the position-velocity slice, but that 
both the model channel map and position- velocity 
morphologies are inconsistent with the data. 

We thus proceed to model NGC 6503 with both 

a thill and a thick layer. As in Schaap. Sancisi. &: Swaterj 
()2000l) . we find that plausible models with co- 
rotating thick and thin layers have similar mor- 
phologies to our one-component models. We 
therefore explore models where the thin and thick 
layers have the same column density profile shape 
and rotation profile shape, but where the thick 
layer has a total column density and rotation ve- 
locity that are lower than that in the thin layer. 

Fig. [TSt l shows one of our better models of 
NGC 6503 with this morphology. The thin layer 
has a FWHM of 0.5 kpc. The thick layer has a 
FWHM of 4.5 kpc, and rotates 25 km slower 
at all R compared to the thin disk. The thin:thick 
layer column density ratio is 4:1. We see that the 
morphology of the position-velocity slice within 
the optical disk approaches that seen in the data. 
However, the model also predicts detectable low- 
velocity gas at larger R in the position-velocity 
slices as well as large-scale, low-intensity emission 
in the channel maps which are not observed in 
the data. Increasing the column density ratio be- 
tween the HI layers alleviates the discrepancies in 
the outer disk and channel maps, but the mod- 
els then fail to reproduce the "beard" across the 
optical disk for plausible differences in velocity be- 
tween the layers. We therefore find that this class 
of models, too, fails to reproduce the position- 
velocity slice and channel map morphology at all 
locations in the HI disk. 

Finally, we invoke 2-component models where 
the thick, low column-density, lagging HI layer is 
truncated beyond the optical radius. One of our 
better models with this morphology is shown in 
Fig. [IS?. The thick layer has a FWHM of 4.5 kpc, 
and rotates 25 km slower at all R compared 
to the thin disk. Within the optical radius {R = 
3.55'), the thin:thick layer column density ratio is 
4:1, and the thick layer has zero column density 
farther out. This model is the only one in FigfTSl 
that reproduces the morphology of the major axis 



position- velocity slice as well as the channel maps. 
While not shown here, this is also the case for the 
advancing half of the disk, as well as for other 
channels and slices in the cube. 

Given the sensitivity of our data, we estimate 
that truncated 2-layer models with combinations 
of column density ratios between 4:1 and 6:1, ve- 
locity lags between 20 km s~^ and 40 km s~^, and 
thick layer FWHM between 3 kpc and 5 kpc could 
reproduce the data. However, the column density 
ratio between the thin and thick disks outside the 
optical disk must be >10:1 in order to escape de- 
tection at the 3-5cr level. We therefore conclude 
that the extra-planar HI layer in NGC 6503 ex- 
tends only out to its optical radius. 

5. Discussion and conclusions 

We have presented aperture synthesis obser- 
vations of the nearby, late-type spiral galaxy 
NGC 6503. Advances in calibration and imag- 
ing software enabled us to produce a wide-field 
continuum image at a dynamic range of 9000:1. 
We detect HI emission at column densities as low 
as 1.8 X 10^^ cm~^, with a spatial resolution of 0.35 
kpc at the galaxy (Fig. [5]) : our images thus have 
7 and 6 tim es higher sensitivity and r esolution 
than those o f ivan Moorsel fc Wellsl (Il985l ). respec- 
tively. The HI cube is also to a great extent free 
of contamination due to missing short spacings. 
In Appendix [Aj we discuss the improved imag- 
ing algorithms available within the AIPS software 
package which have made these images possible. 

Even at the higher sensitivity and resolution of 
our observations relative to previous studies, we 
find that the HI disk in NGC 6503 is remarkably 
regular (Figs. O [5]). Nonetheless, several devia- 
tions from an idealized distribution are apparent: 
the outer portions of many emission "wings" bend 
away from the major axis in the channel maps, and 
the zeroth moment reveals a number of density 
enhancements giving the impression of large-scale 
spiral arms as well as smaller spiral-like structures 
closer to the dynamical center. The HI density 
is reduced at the dynamical center and at several 
other locations in poorly-resolved holes about 0.3 
kpc in diameter. 

We trace the rotation curve of NGC 6503 out 
to 3.7 optical radii (Fig. [TT|l . and find that it is 
flat almost out to the last measured point. Fit- 
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ting rotation models to the first moment image by 
quadrants reveals kinematic features comparable 
in size to the channel width (5.15 km s~^), some 
of which probably reflect spiral density wave dis- 
turbances (Fig. [T^ . This is also evident in our 
CUBIT models of the HI density distribution and 
kinematics (Fig.[T3|): the velocity field is distorted 
from the fit model in the SW portion of NGC 6503, 
partially confirming the presence of the warp in 
the pl ane suggested by IShostak. WilHs. fc Grand 
( 198l[ ). However, a comparable disturbance to the 
kinematics is not observed on the opposite side of 
the plane. In addition, the HI density is greatly 
enhanced compared to, with a significant veloc- 
ity discrepancy from, the fit model along the SE 
major axis. Careful consideration of these asym- 
metries will be taken in forthcoming dynamical 
models. 

Despite the greater HI emission-line sensitivity 
of some recent observations of nearby spiral s com- 
pared to those presented here (e.g.fFraternah et al 



2OO2I : lOosterloo. Fraternah. fc Sancisil 120071 ) . we 
find that the disk of NGC 6503 extends to a 
similar lim iting column den s ity (F ig. [8|). In the 



images of iFraternali et al.l (|2002i ). the disk of 
NGC 2403 truncates at ~ 2 x 10^^ cm'^ (see 
their fig. 6) , and the lowest contour in the 



Oosterloo. Fraternali. fc Sancisil ( 2007 ) map of 
NGC 891 is at 1.0 x lO^^cm'^ (see their fig. 1; 
they do not discuss an actual examination of the 
spectra at such points). We trace the HI disk 
in NGC 6503 out to 1.8 x lO^^cm^^ (22.6 kpc). 
But the absence of absorption towards the quasar 
1748-1-700 yields an upper limit of 5 x 10^'^ cm"^ 
for the column density of cold HI gas along a line 
of sight which should intersect the disk at 29 kpc 
(Fig. [7|). Our results therefore support t he re - 
mark by lOosterloo. Fraternali. fc Sancisi (2007) 
that significant improvements in sensitivity have 
not increased the observed radial size of HI disks. 
Perh aps we are observing the cutoff suggested by 
Silk fc S unvacv (1976) due to ionization by an 
intergalactic UV radiation field. 

While disks do not "grow" radially when ob- 
served at higher sensitivity, nearly all exhibit 
extra-planar HI at lower velocities than those in 
the plane of the disk (see Section[l]for references). 
Our images, when smoothed to achieve greater 
sensitivity, are no exception (Fig. [T5|). We model 
the vertical extent and kinematics of this extra- 



planar HI, and find that the channel maps and 
position- velocity images of NGC 6503 are best de- 
scribed by 2-component models with a thin disk 
and a low column density, lagging thick disk that 
is truncated at the optical radius. We estimate 
that 2-component models with combinations of 
thin:thick layer column density ratios between 4:1 
and 6:1, velocity lags between 20 km s~^ and 40 km 
s~^, and thick layer FHWM between 3 kpc and 
5 kpc could reproduce the data. However, the 
thin:thick column density ratio beyond the opti- 
cal radius must be > 10 : 1 for consistency with 
our images. 

There is therefore a clear correlation between 
the presence of extra-planar gas in NGC 6503 
and its stellar disk. This morphology is hard 
to reconcile with models in which accreting gas 



forms the thick HI disk ( Barnabe et al. 2006t 



Kaufmann et al.ll2006h . NGC 6503 therefore pro 
vides clear evidence for a galactic fountain origin 
for most of the extra-planar gas in nearby galax- 
ies, as sugges ted by the models of NGC 2403 and 
NGC 891 bv IFraternali fc BinnevI (|2008l ). If 10- 
20% of the extra-pla nar gas does stem from accre - 
tion as suggested bv IFraternali fc Binnev (2008), 
it may be possible to detect its interaction with the 
thin disk beyond the optical radius of NGC 6503. 
Considering our limits on the column density ra- 
tio between the thick and thin layers, observations 
with 2-3 times the sensitivity of those presented 
here would be needed to detect this phenomenon. 
Such sensitive observations would also allow for a 
detailed search of High Velocity Cloud analogs in 
a galaxy with a modest star formation rate, which 
may shed some insight into their origin as well. 
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A. Image Deconvolution in IMAGR 

The imaging and deconvolution of most data collected by the VLA are performed with the AIPS task 
IMAGR. Since a number of algorithms within this program have a significant impact on the quality of the 
resulting images, we elaborate on its functionality here. We focus on the multi-scale Clean implementation 
in IMAGR, and illustrate some of its advantages over the standard point-source Clean using the NGC 6503 



data p resented in this paper. A similar implementation of the multi-scale Clean is des cribed by Cornwell 
(|2008l ) and an extensive study of its use in imaging HI in galaxies has been described bv lRich et al.l (|2008l ). 



Their conclusions regarding this method are similar to ours. 

A.l. Cotton-Schwab-Clark Clean 

Clean was proposed by iHogb omi (jl974l ) as a method to determine the sky brightness distribution of 



an object from incomplet e ly sam pled visibilities. A review of Clean and its variants can be found in 
Cornwell. Braun fc Brigg^ ( 1999( ): briefly, a suite of suitable models are convolved by the beam synthe- 



sized from the visibility data and iteratively subtracted from an image until a user-specified criterion is 
reached. The model components are then convolved with a "Clean" representation of the synthesized beam 
(typically a Gaussian) and added back into the residual image to form the sky brightness distribution. 

When the Clean algorithm works solely in the image domain, the convolution of the m odel w i th the 
synthesized beam must be carried out over an area four times larger than that to be Cleaned. IClark (|l980l) 



pointed out that a smaller "beam patch" could be used to find a modest number of model components, after 
which the residual image could be re-computed by subtracting the current source model fully. Thi s image- 
based Clark Clean did the subtraction via Fourier transforms of images. Cotton and Schwab (see SchwabI 



19841 ) were the first to suggest and implement an extension of this method. In the "Cotton-Schwab-Clark" 
version of Clean, after the modest number of model components are found, their visibilities are subtracted 
from the residual visibility data and a new residual image is formed; the process then continues. This allows 
model components to be found over most of the image while avoiding Fourier transform aliasing and other 
inaccuracies in the model subtraction. Cumulative computational error in the imaging is also reduced since 
the images are continually recomputed using the current residual data which gradually approaches the noise 
levels. The Cotton-Schwab-Clark Clean is the default deconvolution algorithm in IMAGR. 

Unconstrained, th e Clean algorithm can move flux density from real objects into erroneous locations 
( Condon et al. The synthesized beam of well-sampled visibility data from the VLA reduces this 



"Clean bias," because of the relatively low sidelobe levels. To reduce it further. Cleans are typically restricted 
interactively ( "boxed" ) to those areas clearly containing real emission. Any emission regions too weak to 
be included will have negligible sidelobes and, thus, do not need to be deconvolved. We compare the flux 
densities measured from images deconvolved with unrestricted and boxed Cleans for the particular case of 
NGC 6503 in Appendix Ol 

A. 2. Wide-field imaging 

When images are constructed with two-dimensional Fourier transforms from interferometers that are not 
co-planar, errors occur which increase as the square of the angular distance from the phase reference point. In 
order to avoid serious distortions in source sha. pes and coordinates, images of large areas must be made up of 
rather smaller facets (jCornwell fc Perlev Il992h . The choice of facet location is significant, since aliasing and 



other numerical errors render the edges and, especially, the corners, of images unreliable. AIPS implements 
an algorithm due to W. D. Cotton in which the facets are placed on circles so that the full area may be 
Cleaned without using the facet corners. 

When imaging a particular facet, IMAGR rotates the projected baseline coordinates and the observed 
phases to the phase reference position at the facet center. An inverse rotation is performed when the model 
components in the facet are subtracted from the current residual data. These rotations are simple matrix 
multiplications and and generally have a negligible computational cost. Additionally, the time required 
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to deconvolve correctly-rendered emission is considerably less than that required to attempt to deconvolve 
distorted renderings. The computed synthesized beam varies slightly from facet to facet. IMAGR differs 
from other faceting implementations in that it selects one facet to be Cleaned in the current cycle, finds new 
model components in that facet, and subtracts them from the current residual data. When it selects the 
next facet to Clean (see Section lA.Sp . it re-images that facet using the latest residual visibility data. In this 
way, sidelobes of a strong source in one facet are therefore removed from the others before the latter are 
Cleaned, which avoids a significant bias in the final image that is otherwise hard to mitigate. 

A.3. Multi-scale Clean in IMAGR 

Clean was originally descr ibed as a method to model a sky brightness distribution as a collection of 



point sources (|Hogbomlll974l ). While a point-source Clean deconvolves poorly-resolved objects accurately 



and efficiently, it requires a large number of compone nts to represent more extended objects and those 



representations can have non-physical attributes (e.g., ICornwelll 119831 ). In particular. Clean has a bias 
against pixels adjacent to one recently selected. To mitigate artifacts in the final image from this bias, model 
objects are restored to the residual image only after convolution with a "Clean beam." Nonetheless, images 
of extended objects deconvolved wi th a point-source Clean tend to show systematic corrugations that can 
be confused for real structure (e.g.. ICornwell. Braun fc Briggj 19991 ) . 



As early as 1974, one of us (EWG, unpublished) began investigating the use of Gaussian source models 
in the Clean algorithm instead of point-source models. Extended sources of a single size cannot accurately 
model point sources, however, and are of little interest in a n astronomical context sin ce most fields contain 
unresolved objects. In order to mitigate these difficulties, Wakker fc Schwa"r3 ( 1988f ) developed an image- 



based algorithm which simu l taneou sly solved for point components and Gaussians of a single pre-determined 



size. 



Holdawav fc Cornwelll (|l999f ) described an extension of this algorithm which would model the image 



with Gaussians of a few selected widths as well as with points. Their algorithm was an image-only Clean 
with "cross beams" for subtracting components of one scale from images made for another scale. At each 
iteration, a component was found in the sky brightness distribution imaged on a particular scale, and then 
subtracted from the distributions imaged on all scales. It was said that this algorithm could run with high 
loop gai and almost no user-set steering parameters. 

The multi-scale C o tton-S chwab-Clark Clean implemented in IMAGR is rather different from the 
iHoldawav fc Cornwelll (Il999l) algorithm , and has now been successfully applied to a number of datasets 



(e.g. lThilker. Braun. fc Waherbos 2002, : IClarke fc Ensslin 2006; .Moniian et al..,2006l : iKent et al.ll200S) . The 



user specifies the number of circular Gaussians to be used as model components, as well as their widths. 
Although it is not required, one of these Gaussians should have a width of zero to model unresolved objects 
and other small-scale structures accurately. IMAGR converts these widths into tapers (Gaussian functions 
which lower the weights at larger projected baseline lengths) for each scale and uses them to make images of 
the beam and the data. These tapered beam images are used to correct the data images to units of Jy/beam 
and to determine the Gaussian Clean beam for each scale. IMAGR then discards the beam images and 
constructs new beam images at each scale representing the convolution of the tapered beam with the Clean 
beam. This process is carried out rather simply by re-making each beam image with a taper \/2 smaller, 
where the taper parameter in AIPS is the baseline length at which the tapering Gaussian function reaches 
30%. 

IMAGR Cleans each resolution for each facet separately, treating M spatial facets times N scales as MN 
total facets. At each major cycle, it selects one of the MN facets to Clean, images that facet with the 
current residual visibility data, finds Clean components in that facet, and subtracts the Fourier transform of 
those components from the residual visibility data. It stops Cleaning a particular facet when the peak flux 
density therein falls below a scale-dependent cutoff, and the cleaning stops altogether when all facets satisfy 
this or another termination criterion. The components on all scales are then restored to the image made at 



the amount of each component subtracted at each iteration, usually 0.1, but values close to 1.0 were said to work well 
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the finest (usually point-source) scale. These components are added back as Gaussians of the size fitted to 
the synthesized beams and scaled by them to produce an image with correct units. 

This algorithm does not work well without some steering. The user may specify a variety of multi-scale 
Clean control parameters in IMAGR, but only two of them have turned out to be important. The first is the 
resolution-dependent flux density cutoff discussed above, passed to IMAGR via the parameter FGAUSS. In 
general, lower resolution images should have higher flux density cutoffs measured in Jy per the beam area 
of that resolution. The second important steering parameter controls the order in which facets are Cleaned. 
When choosing which facet to Clean, IMAGR selects the one with the largest weighted peak flux density 
F ■ 

where F is the peak flux density of the facet, R is the ratio of the area of the Clean beam for the facet being 
Cleaned to the area of the smallest Clean beam, and 6 is a user-specified "bias" parameter that is passed to 
IMAGR via the parameter IMAGRPRM(ll). The general idea is to force IMAGR to select different scales 
with equal probability in the next major cycle. Thus, if the sky distribution is entirely composed of point 
sources, F will be the same at all scales and the user should set h — Q. In the other extreme, if the sky 
distribution is very large and fills the image, then F in Jy/beam will be proportional to the beam area of the 
facet being Cleaned, and the user should set 6=1. In practice we find that 0.2 < b < 0.7 is appropriate for 
data obtained at the VLA, the precise value depending on the structure of the sky brightness distribution. 

Multi-scale Clean tends to be faster than point-source Clean despite the added number of facets. The lower 
resolution facets Clean out a large portion of the total flux density with a modest number of components, 
leaving the high-resolution Clean to image point sources and the edges of objects left behind by the circular 
Gaussian models. Because it makes better use of the short-spacing data, multi-scale Clean also tends to 
reduce the negative "bowls" su rrounding sources observed with missing short-spacing visibility data (e.g., 
Cornwell. Braun fc Briggdll999l ). This is illustrated in the specific case of NGC 6503 below. 



A.4. Example: NGC 6503 

The images of NGC 6503 presented in this paper were produced with IMAGR's multi-scale Clean algorithm 
(see Section[2]). Below, we use these data to illustrate the advantages of the multi-scale Clean over the default 
point-source Clean. 

The continuum in the vicinity of NGC 6503 was imaged and deconvolved with a multi-scale Clean over 
an area 2.5° in diameter divided into 55 facets, using Gaussian source models of FWHM O"(point source), 
36" and 108", flux density cutoffs of 0.21, 0.42 and 1.7 mJy/beam respectively and b = 0.62. The difference 
between the image in Fig. [1] and that obtained with a point-source only Clean for this dataset is shown in 
the left-hand side of Fig. [161 The boxes used in the two Cleans were identical. Because the continuum in 
NGC 6503 is extended (Fig. [2]), the multi-scale Clean managed that area of the image rather better than 
the point-source Clean, and the difference between the two images at this location peaks at ~ 60 /zJy. The 
rest of the continuum image is dominated by emission from unresolved objects (Fig.[T]): for this distribution, 
the benefits of a multi-scale Clean over a point-source Clean are minimal. Nonetheless, the difference image 
in Fig. [16] shows a circular pattern surrounding NGC 6503, with peaks in excess of ~ 20 /iJy. This circular 
pattern is the Fourier transform of the roughly circular hole in the data sampling due to missing spacings 
less than the antenna diameter. Multi-scale Clean appears to have removed a systematic defect from the 
image, reduced the noise by a few percent, and raised the flux density on the extended source by ~ 9%. 

The HI in NGC 6503 was imaged in each channel in a single facet, and deconvolved using a multi-scale 
Clean with 0", 36", 108" and 324" FHWM source models, flux density cutoffs of 0.2, 0.55, 2.6, and 7.8 
mJy/beam respectively and b — 0.62. Because the channels in Fig. [3] contain extended emission, a multi- 
scale Clean performs substantially better than a point-source Clean for this cube. As an example, multi-scale 
Clean found 16% more flux density in 6% more iterations than the point-source Clean at V = 98.3km s~^ 
using identical boxes; see the right-hand side of Fig. [TH] Relative to a point-source Clean, the multi-scale 
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Clean of the HI data cube has ehminated most of the frequency-dependent "bowl" of negative emission 
indicative of missing short spacings in the visibility data. This is illustrated by the primary beam-corrected 
integrated spectra of NGC 6503 in Fig. [171 The dashed line in the left panel shows the flux density of 
NGC 6503 Cleaned and measured in a rectangular region encompassing the zeroth moment image (Fig. [5]). 
It falls systematically below the spectrum measured from the blanked data cube discussed in Section 13.21 
(solid line, reproduced from Fig. llOp because negative emission from a bowl has been included in the sum. 
By contrast, the dotted line shows the same flux density measurement as the solid one, but for a point- 
source Clean data cube of the same visibility data with the same zeroth- moment, channel- independent box 
area. The dotted line falls well below the dashed one because the bowl is much more pronounced in the 
point-source Cleaned data cube than in the multi-scale Cleaned one: 33% less flux density is measured from 
the former in this region. In the right-hand panel of Fig. [ITl the same three integrals are plotted, but all 
three are restricted to the un-blanked regions of the cube that was carefully boxed during, and carefully 
blanked after, Cleaning. The loss of flux density due to the negative bowl has been eliminated for the 
multi-scale Clean (total raised from 185 to 200 Jy km s~^), but a loss of ~ 10% remains for the point-source 
Clean (total raised from 124 to 182 Jy km s~^). The cube which was Cleaned using a single large box area 
encompassing the full HI galaxy, but which was carefully blanked thereafter (solid line) has 3 Jy km s^^ 
more total flux than the carefully Cleaned cube (dashed line), perhaps partly due to Milky Way emission 
and partly to regions of emission overlooked in the careful Cleaning. The general conclusion of this is that 
careful, channel-dependent boxing during the multi-scale imaging in not required, although some restriction 
of the Clean area is always recommended. Careful blanking of the resulting image cube will then eliminate 
noise regions and some remaining negative bowl. Note that, if there is any negative bowl, then the resulting 
total fluxes will be reduced by the fact that signal regions are sitting inside the bowl and so have too low 
a value. Since the negative bowl problem is very much worse with the point-source Clean, the reduction of 
the brightness within signal regions remains significant. 

The disparity between the left- and right-hand sides of Fig. [T7] serves as a reminder that a multi-scale 
Clean — as with any other deconvolution algorithm - does not preclude the need for a careful analysis of 
the data. Nonetheless, IMAGR's multi-scale Clean algorithm outperforms a point-source Clean for extended 
sources like NGC 6503, yielding a truer representation of the sky brightness distribution. 
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Fig. 16. — In gray-scale, the multi-scale Clean image minus point-source Clean image of (left) the continuum 
emission and (right) the HI line emission a.tV = 98.3 km s~^ in NGC 6503. The gray-scale ranges are shown 
in the step wedges, labeled in /xJy. The multi-scale; images in the continuum and the HI are sketched in with 
contours. The Clean beam is illustrated at the lower left corner in both images and is nearly invisible on 
these scales. 
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Fig. 17. — Primary beam-corrected integrated spectrum of NGC 6503 obtained using a number of Cleaning 
and blanking prescriptions. Left: the solid line is the integral over the blanked cube, made from a carefully 
boxed Clean of each channel image (reproduced from Fig. [T0|) . The dashed line is the integrated spectrum 
over a rectangular region slightly larger than the extent of the zeroth moment image constructed with multi- 
scale Clean with little restriction of the Cleaning area. The dotted line is the integrated spectrum constructed 
in the same way but using only point-source Clean components. Right: Same as left, but only un-blanked 
regions in the carefully boxed cube are included in each profile. 
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